Frustrated Quantum Spin Models with Cold Coulomb Crystals 
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We exploit the geometry of a zig-zag cold-ion crystal in a linear trap to propose the quantum simulation of a 
paradigmatic model of long -ranged magnetic frustration. Such a quantum simulation would clarify the complex 
features of a rich phase diagram that presents ferromagnetic, dimerized antiferromagnetic, paramagnetic, and 
floating phases, together with previously unnoticed features that are hard to assess by numerics. We analyze in 
detail its experimental feasibility, and provide supporting numerical evidence on the basis of realistic parameters 
in current ion-trap technology. 



Our understanding of interacting quantum many-body sys- 
tems is usually hindered by their inherent complexity. In or- 
der to unveil their properties, an original approach are the 
so-called quantum simulations (QSs), which exploit a well- 
controlled quantum system to study a complicated interacting 
model [1]. In particular, the QS of magnetic ordering is a 
promising direction of research. Although this effect was al- 
ready identified by the ancient greeks, it still presents some 
puzzles that would benefit from a quantum simulator since 
they cannot be efficiently addressed numerically [2]. Among 
the most prominent platforms for the QS of magnetism, such 
as neutral atoms or photons [3], trapped ions [4] have the ad- 
vantage of near unity fidelity in state preparation, individual 
readout, and a wide tunability of interactions. Here, we in- 
vestigate the QS of frustrated magnetism with trapped ions, a 
phenomenon that yields a fundamental challenge with impor- 
tant connections to high-temperature superconductivity and 
other exotic states of matter, such as spin liquids/glasses [2]. 

Even if the typical distances between trapped ions make 
a direct dipole-dipole coupling negligible; the long-range 
Coulomb interaction yields an indirect mechanism whereby 
the vibrational phonons, coupled to the spins by lasers, me- 
diate a spin coupling [5]. In the seminal work [6], frustra- 
tion in a 3 -ion chain has been achieved by fine tuning the 
laser frequencies within the collective vibrational band. How- 
ever, scaling this scheme to larger ion chains may find a fun- 
damental limit: the vibrational band becomes increasingly 
dense, and the lasers approach a certain resonance. An alter- 
native scheme would use micro-fabricated surface traps [7], 
but the current ion-heating problems must be overcome to 
obtain strong-enough interactions. One may instead rely on 
large trapping non-linearities [8], yet the requirements remain 
to be met in experiments. Here, we present a different method 
amenable of being scaled to larger ion numbers that relies on 
state-of-the-art ion traps and spin-dependent light forces [9] . 

Trapped ions crystallize in self-organized Coulomb crys- 
tals at low temperatures. By tuning the trapping frequen- 
cies [10], a linear string of ions undergoes a structural tran- 
sition to a two-leg zig-zag ladder [Fig. 1(a)]. Whereas the 
linear configuration has been the usual scenario for the QS of 
magnetism [4], the capabilities of the zig-zag ladder have re- 
mained largely unexplored. In this Letter, we fill in such gap. 
We show that by exploiting the commensurability of the inter- 




Figure 1. Zig-zag ladder: (a) Equilibrium positions in units of l z for 

N = 30 ions, COy > co x = 6.lco z , as given by d(V t rap +Vcoul)/^0'a = 
0. (b) Scheme of the stimulated Raman transition. A pair of 
lasers with frequencies &>i,C02, and wavevectors ki,k2, are far- 
detuned from the transitions and lead to a two-photon transition 
ft>L = (0\ — o>2. (c) Transverse, Q.„, and planar, Q!^, normal-mode 
frequencies. In the inset, the amplitudes of the highest-frequency 
modes. We also show the frequency &>l of a laser beam blue-detuned 
from the vibrational band leading to the transverse spin-dependent 
force, (d) Scheme of the spin interactions, where J\ < {J 2 > 0) 
is a ferromagnetic (antiferromagnetic) coupling. Incompatible frus- 
tration of two plaquettes due to the dipolar range, (e) Coupling be- 
tween the central ion 7*0 = N/2 and its neighbors, as obtained for 

* l ,a=*/(4-4 +1 ). 



leg distance of the ladder with the wavelength of the spin- 
dependent light force, it is possible to control both the relative 
strength and the sign of the spin interactions in a way that is 
almost independent of the crystal size, and can be thus scaled 
to larger ion crystals. We describe how to reach regimes of 
competing ferromagnetic (F) and antiferromagnetic (AF) cou- 
plings, which are incompatible with the crystal geometry so 
that the groundstate cannot minimize all interactions simul- 
taneously. Therefore, trapped ions can be used as quantum 
simulators of the two main sources of frustration, namely, ge- 
ometric frustration and competing interactions. 

Let us summarize our main results: (i) The sign of the spin- 
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spin interactions between distant ions depends on wether they 
lie along the same or different rungs of the ladder. ( ii) The 
relative strength of these two interactions can be controlled by 
tuning the ratio of the inter-leg distance to the wavelength of 
the light forces. (Hi) These ingredients yield a fully-tunable 
and experimentally-feasible quantum simulator of the /1-/2 
quantum I sing model (J1-J2-QIM), a cornerstone of magnetic 
frustration [11]. (iv) The long-range character of the spin cou- 
plings introduce incompatible sources of frustration leading to 
a richer phase diagram with a novel conjectured phase. 

Model- We consider N ions of mass m and charge e, con- 
fined in a linear trap with frequencies {co a }a=x,y,z [12]. The 
geometry of the Coulomb crystal r®, determined by the bal- 
ance of the trapping forces and the Coulomb repulsion, cor- 
responds to a zig-zag ladder in the x-z plane for co z < (Q x <C 
COy [10]. The vibrations around the equilibrium rj = r*- + Ary 
are coupled through the Coulomb interaction, which yields 
model of coupled oscillators in the harmonic approximation 
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where are the vibrational couplings [13]. We identify 
two types of collective excitations, namely transverse and pla- 
nar modes that describe the motion perpendicular and parallel 
to the crystal in terms of the phonon operators a nJ b n . The vi- 
bration amplitudes jtfij- n ,jtfi"- n , and frequencies are 
obtained by solving the quadratic problem (1), leading to 
#ph = L^=i ^n a n a n + £n=i ^nHi^n w i m me phonon branches 
displayed in Fig. 1(c). As shown in this work, when the pla- 
nar bandwidth is small enough, a blue-detuned laser beam will 
only excite the transverse phonons even if it also has a compo- 
nent along the crystal [Fig. 1(a)]. This is a crucial point of our 
scheme. The transverse phonons, which are best suited to act 
as carriers of the magnetic interaction, will be responsible of 
mediating the spin- spin coupling. The component of the laser 
along the zig-zag plane will not excite the planar phonons, 
but rather open the possibility of tailoring the ferromagnetic- 
antiferromagnetic sign of the spin coupling depending on the 
ratio of the inter-leg ion distance with the wavelength of light. 

Effective frustrated spin models - We consider a pair of 
laser beams that induce a stimulated Raman transition be- 
tween two internal states |t/)?|l/) [Fig- 1(b)]. When 
the difference between the laser frequencies matches the 
transverse vibrational frequency, we get a Stark shift = 

^ly^^^'^r^+Ec., where oj is a Pauli matrix, 
is the two-photon Rabi frequency, and kL lies in the x-y 
plane. After introducing the phonons, we make an expansion 
for small Lamb-Dicke parameters 7] w _l = &L,y/ V^ m ^n ^ 1 
(similarly for the planar modes). By tuning the laser frequency 
above the transverse phonon branch, we can neglect all non- 
resonant terms but the transverse spin-dependent dipole force 
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where S n ± = £1^ ~ and we work in a frame rotating with 
the phonons. For this approximation to hold, the following 
conditions must be satisfied: <C col, tj wm Ql <^ \&n ~ ^1- 
The first condition is necessary to suppress the carrier terms, 
and the second one to avoid the coupling of the spins to the 
planar modes, and are both satisfied for the realistic parame- 
ters below. Additionally, by aligning the lasers nearly parallel 
to the v-axis, we favor the fulfillment of the last condition. 

The conditional force (2), which couples the spins to the 
transverse phonons, induces a spin interaction due to the vir- 
tual exchange of phonons. The exact expression is obtained 
by a Lang-Firsov-type transformation [14], [5] that decouples 
spins and phonons e _iS (// p h +H^)g s w // e ff + if p h, where S = 
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/ j n a z -a\ — H.c. generates a spin-dependent 
displacement. This leads to an effective Ising model 
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where the transverse field, h, follows from a microwave res- 
onant with the transition, and Gj is a Pauli matrix. In this 
work, we are interested in the anisotropy of J^f oc cos(kL • r^), 
which has been overlooked so far since it vanishes for linear 
ion chains. Remarkably, it can lead to magnetic frustration in 
planar crystals. In the limit of tight transverse confinement, 
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where T] a = &L,a/ \j2mQd a , and l z = (e 2 /mco 2 ) 1 ^ 3 gives a unit 
for the distances. A crucial observation for the QS of frus- 
tration is that the sign of the couplings may be different for 
ions that belong to the same rung of the ladder, or to distinct 
ones. More precisely, if the condition k^J z = nn/ (x® — ^+1) 
is fulfilled for a pair of ions in the bulk of the crystal such that 
n is an integer, we obtain a ferromagnetic (antiferromagnetic) 
coupling between pairs of ions in different (same) rungs. We 
corroborate this claim by obtaining the exact couplings from 
the numerical solution of Eq. (3), which yields the aforemen- 
tioned sign alternation [Fig. 1(e)] in agreement with the ap- 
proximate description (4). Let us highlight the importance of 
this result: the groundstate shall not be able to minimize such 
spin-spin couplings simultaneously, and will thus yield a frus- 
trated magnet [Fig. 1(d)]. According to the Tolouse- Villain 
criterion [15], our spin model (3) is frustrated since there is an 
odd number of AF couplings per unit cell. In the following, 
we support this analytical derivation with numerical simula- 
tions on the basis of the realistic experimental settings and 
parameters that are described now. 

Experimental setup- To realize the spin model (3), we rely 
on trapped-ion quantum-processing technologies developed to 
a high degree of perfection [16]. We emphasize that all pro- 
cedures and experimental parameters are well in reach with 
current experimental standards. Cold-ion crystals are kept in 
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Figure 2. QS of the frustrated quantum Ising model: (a) Dynamics of (Oj(t)) for an initial spin excitation |+)i in the allow ed-hopping 
regime. Comparison of the analytical description , and the numerical study of the complete Hamiltonian -\-H^ where the phonon Hilbert 
space is truncated to one excitation, (b) In the inhibited-hopping regime, there is a complete agreement between both descriptions, where the 
spin excitation cannot occupy site 2. (c) Schematic phase diagram of the J1-J2-QIM with ferromagnetic (F), dimerized antiferromagnetic (AF), 
paramagnetic (P), and floating (FP) phases, (d) Conjectured phase diagram with dipolar interactions. Also highlighted in grey the regions 
studied in (e)-(f). (e) q = n/2 component of the structure factor evidencing the AF-FP-P transitions, (f) Different magnetic modulations at 
h = showing the splitting of the critical points with a different order in between, (structure-factor components are normalized with respect to 
the maximum). 



linear Paul traps with a combination of static and dynamic 
electric fields [12]. The specific case in Fig. 1(a) requires a 
significant anisotropy of the radial confinement COy, which 
can be achieved by deforming the symmetric four-rod trap 
electrode configuration. Additionally, the anisotropy can be 
increased and tuned to the desired value by applying a bias 
voltage to the DC electrodes. In this way, the desired zig-zag 
crystal is trapped with MHz trap frequencies, and fulfilling 
the condition co z < (0 X <C COy. The spins in the generic level 
scheme in Fig. 1(b) might correspond to the hyperfine states of 
9 Be + , 43 Ca + ions, or to the Zeeman sub-levels of the ground- 
state of 40 Ca + . These levels are coupled with Raman laser 
beams that mediate the optical spin-dependent light forces [9] . 

The experimental sequence consists of a number of steps 
starting with the ( i) initialization of the system. Here, the tem- 
perature of all modes is reduced in a Doppler cooling step, 
followed by a sub-Doppler cooling of the transverse mode. 
The relatively high frequency (Qy/2n — 20MHz ensures good 
cooling results. The spins are initialized by optical pumping, 
followed by a n/2 pulse on a resonant Raman transition (not 
shown in Fig. 1(b)) such that a spin superposition is reached. 
(ii) The spin-spin interaction is generated by application of 
the Raman light fields. To suppress spontaneous processes, a 
large detuning A > 10GHz is chosen, and light beams with an 
optical power of about lOmW are focused into a waist of 30 
to 60/im [17] suffice to generate / e ff ~ 1kHz. After an inter- 
action time in the millisecond range, ( Hi) the simulation result 



is readout from the fluorescence imaged on the chip of a CCD 
camera. For this, we illuminate the entire crystal with light 
resonant to a dipole-allowed transition such that the 1 1) state 
emits fluorescence while | \) does not. For the demonstration 
of frustration, the requirements are much released in compar- 
ison to those of full quantum- state tomography [18]. 

Numerical validation - We address the simplest scenario 
where frustration occurs: a 3-ion zig-zag with trap frequen- 
cies C0y/ca z = 20, co x /co z = 1.43, co z /2n w 1MHz, leading to 
r°j/l z e {(-0.22,0, -0.92), (0.44,0,0), (-0.22,0,0.92)}. We 
test the crucial assumption underlying the derivation of the 
frustrated Hamiltonian (3), namely the possibility to neglect 
the pushing force on the planar modes. Even if the dipole 
force only acts in the x-y plane, the z-motion gets coupled 
through the Coulomb interaction. Hence, we must solve nu- 
merically the dynamics of the complete Hamiltonian // p h + 
and compare it to the effective analytical description H e ff 
in (3), where we set h = 0, \Yo) = |+)i ® |— )2® | — >3, with 
|±) = (| t) =t I i))/V2, and consider a ground-state cooled 
crystal for simplicity. We study two regimes: (i) Allowed 
hopping: We set T] } =0.1 = 10i] x , (Dl 1.2c0y, = 0A\(Oy- 
(Oi\/r\y, and study the evolution of the whole spin-phonon sys- 
tem for &l, x = 0- Since there is no modulation of the coupling 
sign, the initial spin excitation located at site 1 hops around the 
triangular plaquette as a consequence of the Ising coupling. In 
Fig. 2(a), we show the agreement between the analytical and 
the numerical descriptions, (ii) Inhibited hopping: We use the 
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same parameters, but set &l, x = 7i/2l z (x2 — x®). It follows that 
the Ising couplings between 1-2 and 2-3 vanish jf£ = J%3 = 0> 
and thus the excitation only tunnels along 1 o 3, as confirmed 
by both analytical and numerical descriptions in Fig. 2(b). 

Due to the complexity of the spin-phonon system, the nu- 
merics for ground-state cooled crystals already exhaust the 
capabilities of classical computers. However, by studying 
the corresponding Heisenberg equations, one finds (<Jj(t)) = 
(<7j(?)) e ff(l — e), where the finite temperature leads to £ oc 

LnN^QlnL^ji^ji/ 8 ^ with ^ = (4^)th- Therefore, 
there exists a trade-off between strong spin-spin interactions 
and small residual errors, namely, the larger the laser Rabi 
frequency is, the stronger the spin interactions become but 
also the bigger the residual error gets. We have found a com- 
promise for the two, such that for tight confinement and the 
above parameters, this error is roughly £ w 10 _2 /V, and thus 
ground-state cooling is not required. Another possible source 
of error is the micromotion [12], which may cause heating, 
or modify the spin-dependent force by additional sidebands. 
The heating is negligible far from the structural transition 
where non-linearities are small. Additionally, the effect of 
the micromotion sidebands can be neglected if (Oq ^> £1^ , and 
(Ol ~ C0y <C ^rf , where £l r f is the r.f. trap frequency. 

Frustrated quantum Ising model - To describe the essence 
of our quantum simulator, we review the main features of 
the next- to-nearest neighbor quantum Ising model {J\-Ji- 

QIM), H = JiLjOjOj+ 1 +J2LjOjOj+ 2 ~ h ^j a h in the 
frustrated regime sign(/2) = — sign(/i) = +1. For both 
J2,h — ^ 0, one finds a ferromagnetic two-fold ground- state 
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The F-AF phase transition is known to be of first order with 
an exponentially degenerate ground-state at |/i | = 2/2, which 
is a hallmark of frustration. Making h^> |/i|,/2> the ground- 
state becomes a paramagnet with the spins aligned with the 
transverse field |P) = l^ - ^--^ - ^). The transitions be- 
tween the ordered phases mentioned above and the disordered 
paramagnet yields a rich phase diagram [Fig. 2(c)]. Indeed, 
both numerical [19] and analytical [20] methods hint at the 
existence of a new phase separating the AF-P transition. 
This critical phase, commonly denoted as the floating phase 
(FP), exhibits quasi-long range order in combination with 
an incommensurate periodic modulation of the spin-spin 
correlators that depends on the couplings. We stress that the 
existence of the FP has raised some controversy [21], which 
makes a QS of the utmost interest to settle the problem. 

Our QS naturally includes the effect of longer-range inter- 
actions, which may be expected to play an important role in 
the presence of frustration. Using a Lanczos algorithm opti- 
mized for our dipolar model (3), we have diagonalized sys- 
tems of up to 24 spins, and considered a homogeneous lattice 
spacing with o$ +1 = erf in order to capture the bulk prop- 
erties of the ion crystal. To identify the different phases, 



we have computed the magnetic structure factor S zz (q) = 
Lj,A:( <7 / |) el ^^ _ ^ w i m Q £ [0,27r), which is an order param- 
eter accessible in trapped-ion experiments (see below). In 
Fig. 2(d), we summarize our main results, according to which 
the ordered phases of the next-to-nearest neighbor model sur- 
vive to the longer-range interactions. This is evidenced in 
Fig. 2(e), where we show that the S zz (n/2) component of the 
structure factor captures the dimerized-antiferromagnetic and 
the floating phases. As one crosses the floating phase and en- 
ters in the paramagnetic region, the two different types of de- 
cay are consistent with the AF-FP Pokrovsky-Talapov and the 
continuous FP-P Kosterlitz-Thouless phase transitions. 

Remarkably, our results show that the longer dipolar range 
is crucial in the strongly-frustrated regime \ J\\~ 2Ji- We find 
that the additional couplings split the F-AF critical point, and 
moreover we conjecture that they give rise to a new phase 
[Fig. 2(d)]. The different ordering in this central region is a 
consequence of the competing frustration mechanisms intro- 
duced by the long range [Fig. 1(d)]. By including consecutive 
terms of the dipolar tail, we have checked that the splitting of 
the critical points is caused by the third- and fourth-nearest 
neighbors couplings. By increasing these couplings at the 
expense of the nearest- and next-to-nearest-neighbor interac- 
tions, the size of this central region grows. Besides, its robust- 
ness is not modified by including further longer-range terms, 
such as fifth- and sixth-neighbor interactions. In Fig. 2(e), one 
observes that this new central region supports a different or- 
dering, as evidenced by the existence of a modulation in the 
structure factor with momentum q = 71 / 4 not present in the F, 
and AF phases. We have checked the decay of this component 
as one moves into the P phase, which supports this conjecture. 

Efficient detection methods for the phase diagram .- Since 
full quantum state tomography becomes inefficient for many- 
body systems, we focus on the measurement of the order pa- 
rameters above. One of the advantages of trapped ions with 
respect to other platforms is their ability to perform high- 
accuracy measurements at the single-particle level [12]. From 
the spatially resolved state-dependent fluorescence, one can 
infer local observables mj :Z = (Oj) = ^(PJ — 1), or two-body 

correlators Q,(r°, r°) = <lrfof> = \ [1 - 2(Pj + P e k ) + 
An alternative is to measure global properties of the light 
emitted by the ion ensemble, which carries information about 
the spin correlations [22]. The fluorescence spectrum along 
a detection direction, r, is related to the structure factor 
O) oc I^e^'^-^^l + of)(l + of)), where k = co/c is 
the wavevector of the emitted light. Even if the ion spacing is 
larger than the optical wavelength, one may compensate it by 
placing the detector orthogonal to the crystal plane, and us- 
ing a good angular resolution. Note that such structure factors 
may be used to obtain a lower bound on entanglement [23]. 

Conclusions and outlook- We have demonstrated theoret- 
ically that cold Coulomb crystals can be used as quantum 
simulators of paradigmatic, yet controversial, models of frus- 
trated quantum spins. We presented a method to control the 
F-AF nature of an Ising spin-spin interaction between the ions 
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trapped in a zig-zag structure, leading to a quantum simula- 
tor of the J\-J2 quantum Ising model with dipolar range. We 
have discussed the experimental feasibility of the proposed 
work and corroborated it by numerics based on experimen- 
tally realistic parameters. Analogous ideas can be transferred 
to different geometries in self-assembled Coulomb crystals, or 
surface microtraps, to realize other frustrated models. These 
ideas might also be relevant to atoms confined in multimode 
optical cavities 
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